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ABSTRACT 



Aims. We aim to derive a stability condition for non-isothermal Bonnor-Ebert spheres and compare the physical properties of critical 
non-isothermal and isothermal gas spheres. These configurations can serve as models for prestellar cores before gravitational collapse. 
Methods. A stability condition for non-isothermal spheres is derived by constructing an expression for the derivative of boundary 
pressure with respect to core volume. The temperature distribution is determined by means of radiative transfer calculations. Based 
on the stability analysis, we derive the physical parameters of critical cores for the mass range 0. 1 - 5.0 M Q . In addition, the properties 
of roughly Jupiter-mass cores are briefly examined. 

Results. At the low-mass end the critical non-isothermal sphere has lower central density and a slightly larger physical radius than the 
corresponding isothermal sphere (i.e. one with the same mass and average temperature). The temperature decrease towards the core 
centre becomes steeper towards smaller masses as the central density becomes higher. The slope depends on the adopted dust model. 
We find that the critical dimensionless radius increases above the isothermal value £o = 6.45 for very low-mass cores (< 0.2M o ). 
However, in the mass-range studied here the changes are within 5% from the isothermal value. 

Conclusions. The density structures of non-isothermal and isothermal Bonnor-Ebert spheres for a given mass are fairly similar. 
However, the present models predict clear differences in the average temperatures for the same physical radius. Especially for low- 
mass cores, the temperature gradient probably has implications on the chemistry and the observed line emission. We also find that 
hydrostatic Jupiter-mass cores with radii less than 100 AU would have very high boundary pressures compared with typical pressures 
in the interstellar space. 

Key words. ISM: clouds - Radiative transfer 



1. Introduction 

Physical conditions prevailing in prestellar cores, i.e. gravita- 
tionally bound concentrations of molecular gas, constitute the 
initial set-up of star formation. Prestellar cores are generally 
not spherical in sh ape, but it has been found that a Bonnor- 
Ebert (BE) sphere (lBonnoJll956t lEberdfl955l) . i.e. an isother- 
mal gas sphere in hydrostatic equilibriu m, can succesfully ap- 
proximate several prestellar cores (e .g. iBacmann et al.l l2000l 
lAlves et al.ll200lUKandori et al.ll2005t) . While the density struc- 
ture predicted by a BE model can give a good fit to observations, 
there may be problems with e.g. the p ressure required for the 
core to be stable (see e.g. the review bv lBergin & f^fallal[2007L 
and the references therein). Furthermore, the isothermal nature 
of the observed cores cannot be ascertained, and it is possible 
that at least some of the "isothermal" cores possess a signif- 
icant temperature g r adient due to attenuation of s tarlight (e.g. 
Zucconi et alJ 200 lj_ Ward-Thompson et al.l [2002t iPagani etal] 



20047 ICrapsi et al. I I2007I) . A hydrostatic core with the tempera- 



ture decreasing inwards is one possible configuration preceding 
the gravitational collapse, and can thus serve as the starting point 
for dynamical models. 

In light of the above, it is of interest to study what kind 
of constraints could be imposed on the stability of a non- 
isothermal version of the Bonnor-Ebert sphere (hereafter BES). 
The non-isothermal ("modified") Bonnor-Ebert s phere (here- 
after MBES) has been di s cussed previously by e.g. lEvans et al. l 
(|200lh . iGalli etalJ (I2002I) . iKeto & Field! (l2005l) and lSipila et all 
d2010l) . In this paper, we derive an explicit sta bility condition 
for the MBES, using a method similar to that of iBonnorld 19561) 



in his derivation of the isothermal stability condition. The non- 
isothermal stability condition will allow us to calculate the critial 
radius for the MBES, and this is turn yields information on e.g. 
the possible values of the density contrast and boundary pressure 
for critical modified BE spheres. 

The paper is structured as follows. In Sect.|2] we discuss the 
properties of the MBES and present a method of obtaining the 
non-isothermal solution (£, iff). We also present an explicit sta- 
bility condition and discuss how the critical radius of an MBES 
of given mass can be found. In Sect.[3] we present the results 
of our analysis, i.e. the behavior of the physical characteristics 
(radius; density and temperature structure) of critical cores as a 
function of core mass. In Sect.|4]we discuss the implications of 
the results presented in Sect. [3] In Sect.|5] we present our conclu- 
sions. AppendixlAlincludes detailed derivations of some expres- 
sions presented in Section|2] and AppendixlBl discusses thermo- 
dynamics of the MBES. 



2. The critical radii of non-isothermal Bonnor-Ebert 
spheres 

2. 1 . The modified Bonnor-Ebert sphere 

The Bonnor-Ebert sphere dBonnorf Il956t lEbertl Il955h is an 
isothermal hydrostatic gas sphere in hydrostatic equilibrium. By 
assuming the ideal gas equation of sta te and using the hydrostatic 
equilibrium condition. lBonno n (ll956h obtained an expre ssion for 
the d ensity distribution of such a sphere (eq. (2.3) in iBonnod 
I1956I) . The density distribution equation can be solved by ex- 
pressing the density and radius in terms of non-dimensional 
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variables £ and if/ (eq. (2.6) in lBonnorlll956 ) ; this leads to the 2.2. Solving the non-isothermal Lane-Emden equation 



Lane-Emden equation (eq. (2.8) in iBonnoi Il956l) . which re- 
duces to a pair of first-order differential equations. The solution 
to the Lane-Emden equation then yields the physical parameters 
of the sphere by back-substitution. 

By considering a radial dependence of the temperature, the 
above analysis can be generalized to non-is othermal spheres; 
this leads to a modi fied Bonnor-Ebert sphere dEvans et alj200U 
iKeto & Fieldll2005h . In this case, the density distribution equa- 
tion takes the form 



— — I 

r 2 dr 



dp dT 

T h p 

dr dr 



AnGmp 



(1) 



In ISipila et"aL we discussed a modified version of 

Bonnor's substitutions (eq. (2.6) in lBonno rl 19561) taking into ac- 
count the radial dependence of the temperature: 



Because of its temperature dependence, eq. (0 cannot be read- 
ily solved. One thus needs to supply a temperature profile. In 
order to simulate the physical conditions inside prestellar cores, 
we have chosen to obtain the temperature pr ofile using radia- 
tive transfer modell ing of the dust component (lJuvela & Padoanl 
12003b lJuvelall2005h . We assume in the following that T ias = 
Tdusti this assumption is discussed below. 

Before a temperature profile can be calculated, however, a 
first approximation density profile for the model core is needed. 
A natural starting point for the calculations is the BES, because 
the solution (£, iff) to the Lane-Emden equation is known. 

We first fix the mass of the core. The mass of an isothermal 
core can be written, using the variables introduced previously, as 



M iso = 4n/3^A- 1/2 fo M 



diff\ 



(7) 



p= -e 

T 



R=p m r m^ 



(2) 



(3) 



where r = T/T c , T = T(r) = 7\£), T c is the temperature at 
the core centre, A is the mass density at the core centre and 
f3 = kT c /4-nGm, where k is the Boltzmann constant and m is 
the average molecular mass of the gas (here we assume that the 
gas consists of IT and He, so that m = 2.33 amu). We note that 
our su bstitutions are slightly different than those o flKeto & Field! 
(120051) . who used the temperature at the core edge, r out , as the 
reference temperature. We have chosen the central temperature 
T c as the reference temperature, because it depends less on the 
assumptions of the conditions outside the core. Now, for exam- 
ple, pressure and volume can be expressed in terms of £ and iff 
as 



m 



V = = ^3/2^-3/2^3 



(4) 



(5) 



where we used equations (0 and 01 and assumed the ideal gas 
equation of state p = pkT/m. We note that the dimension in the 
above equations is carried by /? and A (and of course by G), as £, 
and iff are dimensionless. 

Inserting equations (0 and 0) into eq. (0, one obtains the 
non-isothermal ("modified") Lane-Emden equation 



e-2 



d£ 



dip 



1 



(6) 



Imposing the boundary conditions if/ — 0, dif/fdt; — 0, r = 1 and 
drfd^ — at the center one can (numerically) integrate eq. (0 
and thus obtain a density profile for the MBES, although the 
temperature profile must be calculated independently using, e.g., 
radiative transfer means. We discuss this point in further detail 
in the following Section. 

Finally, we note that by choosing T c — T — constant (t = 1), 
one readily recovers the isothermal versions of equations ([1]) to 
(0. 



where j3\ so = kT/4nGm. In the above and from this point on, 
the subscript "out" refers to a value taken at the edge of the 
core. Next, a value for the non-dimensional radius f out has to be 
chosen. The choice of £ out determines the stability of the BES. 
The critical value is £o ~ 6.45; cores with £ out < £o are sta- 
ble, w hereas cores with ^ out > £o are unstable (see e.g. lBoniiorl 
1 1 956b . Also, small values of £ out represent cores with a low cen- 
tral density and rather shallow density gradients, whereas cores 
with large £ out are centrally dense and present a large density 
gradient. 

Having fixed M and f out , we proceed to solve the central 
density A using eq. (0. A density profile is then constructed us- 
ing p = /lexp(-i/'), the isothermal counterpart of eq. (0. A dust 
temperature profile for the core can now be calculated. In the 
temperature calculation, we assume that the core is embedded in 
a larger molecular cloud; we set Ay = 10 at the edge of the core. 
The spectr um of the una ttenuated interstellar radiation field is 
taken from iBlackl (Il994l) . The final calculated temperature de- 
pends on the properties of the dust component. We carry out 
separate temperature calculations assum ing two different types 
of gra ins, using grain opacity data fromlOssenkopf & Henningl 
dl994L henceforth OH94) and from lLi & Draind OOOU hence- 
forth LD01); th e LD01 extinction curve was slightly modified 
as described in Sipila et al ] (I2010h . In the latter case, we con- 
sider both silicate and graphite grains, with grain material den- 
sities 3.5 g cirT 3 and 2.5 g cirT 3 , respectively. From now on, we 
will call core models with OH94 grains "Type 1" models, and 
core models with LD01 grains "Type 2" models. The OH94 and 
LD01 grain models were chosen for this study because these 
models describe grains in the dense medium inside prestellar 
cores, where grain properties and possibly their size distribu- 
tion are thought to d i ffer from those in the diffuse ISM ( e.g. 
ISteinacker et alJl2010l:lPagani et alJl2010l:lJuvela et ap201 ll) . 

In all calculations in this paper, we assume that r eas = 7d us t. 
This assumption should hold well for low mass cores and in the 
central parts of more massive cores (M ~ 4 - 5 M ; hereafter 
these will be referred to as "high mass cores"), but generally in 
the outer parts of high m ass cores r gas + 7d us t (e.g. lGalli et al.l 
120021 iKeto&F ield 2005). However, high mass cores should be 
approximately isot hermal due to low average density (see e.g. 
IKeto & Fieldll2005l) . and hence we can (qualitatively) expect the 
high mass MBES to behave like the BES. 

After the temperature calculation, the modified Lane-Emden 
equation (eq.[6j is solved using the determined temperature pro- 
file. This yields the function iff as a function of The function 
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if/ is not yet self-consistent, however, since the solution was ob- 
tained starting from the isothermal core. We proceed by updat- 
ing the central density of the core using the new function if/. The 
mass of a non-isothermal core can be expressed as 

M noniso =4^ so A- 1 / 2 4 t r out |^J , (8) 

where now /? n oniso = kT c j\itGm. The mass of the core and the 
non-dimensional radius £ out are kept constant during the itera- 
tion, so that the above equation yields a new estimate of the cen- 
tral density A. A new density profile is then calculated according 
to eq. @. A new temperature profile is resolved yielding a new 
solution to eq. (|6j. The iteration is continued until the function iff 
converges; this happens typically after 3-4 iterations. When the 
iteration is complete, the physical parameters of the core can be 
derived using eqs. (f2| to (|5]l- 

2.3. The critical radius of the MBES 

As discussed in Sect. 12.21 the initial choice of £ out determines 
the density structure of the isothermal core, which one uses as a 
starting point for constructing the non-isothermal core. Because 
in the isothermal case the Lane-Emden equation does not de- 
pend on temperature, the solution function ifi is always the same 
regardless of the choice of £ out ; the solution is, in this sense, 
universal. In the non-isothermal case the situation is different - 
each value of £ out represents a unique core configuration with a 
unique temperature structure. This means in particular that the 
solution function if/ depends on £ out . 

The stability condition of both the BES and the MBES can 
be found by constructing the derivative of the boundary pressure 
in terms of core volume. For a stable core, we expect pressure to 
increase in a contraction of the core, corresponding to a negative 
sign of the derivative. The critical point is the smallest value of 
£out for whic h the pressure d erivative turns positive (6p/6V > 0). 
As shown by IBonnorl (I 1 9561) . all values of £ out beyond this point 
represent unstable cores. We thus need to solve eq. © for mul- 
tiple values of £ out , and look for the value of £ ollt for which the 
pressure derivative changes sign. In practice, when the iteration 
is complete, we extract the values of the parameters /3, A and t/Q 
corresponding to each £ out and combine them to form "global" 
functions /3(^ out ), A(i; oal ) and iA(^ ut)- These will be needed to cal- 
culate the pressure derivative. 

The pressure derivative can be constructed by usi ng vari- 
ationa l calculus; our method is analogous to that of IBonnorl 
dl956l) . with the exception that we have extended the discussion 
to non-isothermal models. In practice, the extension is carried 
out by considering variations of three parameters (J3, A and £), 
instead of just the two latter ones needed in the isothermal anal- 
ysis. In the case of the MBES, the pressure derivative takes the 
form (see also AppendixlAb 

W 3V [/3^6f3-A^6A + 2^ QUt \ ' 

We insert into the above equation the global functions f3(£ out ), 
M^om) and i^(^out) and look for the value of £ out for which 
the pressure derivative changes sign; this value is the non- 
dimensional critical radius £i . We have done this for a range of 
core masses; the results of the analysis are presented in Sect.[3] 

1 We extract the value of if/ at the edge of the core, that is at <f oul , so 
that the mass of the core remains fixed. 
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Fig. 1. The H2 number density (solid lines, left y-axis) of critical 
0.25 M Q Type 1 (OH94, blue) and Type 2 (LD01, red) MBESs 
as a function of radial distance from core centre. Also plotted is 
the density profile of the critical 0.25 M Q BES corresponding to 
the mean temperature of the Type 1 MBES (green solid line; see 
text). The temperature profiles of the respective cores are plotted 
as dashed lines (right y-axis). 

The core mass and hence the total number of particles is con- 
served in a contraction of the core, i.e. 6N = 0. This leads to a 
condition between 6/3, 6 A and <5£ out (see Appendix|A]i, which can 
be used to eliminate for example 6/3 in eq. ©■ If this is done, 
it is straightforward to verify that in the isothermal limit, eq. © 
reduc es to the isothermal pressure derivative (eq. 2.16 in lBonnorl 
fl956h . 

3. Results 

We have derived the physical parameters (radius, density, tem- 
perature) of critical MBESs for a range of core masses using the 
iterative approach described in Sect. [2] In this Section we first 
compare the radial density distributions within a critical MBES 
and a critical BES, and thereafter describe how the "global" 
properties, i.e. the central density and the temperature profile of 
an MBES depend on the core mass. 

3.1. Non-isothermal vs. isothermal: internal density structure 

The density distributions of an MBES and a BES for a given 
mass are compared in Fig.Q] Here we plot the gas number den- 
sity (n(H2), solid lines, left y-axis) as a function of radial dis- 
tance from the core centre for a critical MBES of mass 0.25 
M Q . The Type 1 core is plotted in blue and the type 2 in red. 
Also plotted (green solid line) is the density profile of a critical 
BES of the same mass and a temperature which is equal to the 
average temperature of the Type 1 core (cf. eq. IB.4l ). The tem- 
perature profiles are plotted with dashed lines keeping the same 
color codes as for the densities (the scale is on the right). 

On can see in the Figure that the higher central temperature 
of a BES as compared with an MBES permits also a higher cen- 
tral density, and the BES is more compact. A similar effect is 
seen between Type 1 and Type 2 cores. The Type 1 core is 
slightly warmer and more centrally concentrated than the Type 



4 



O. Sipila et al.: On the stability of non-isothermal Bonnor-Ebert spheres 




10 3 r 



itr 1 10° 10 1 

Core mass [M sm ] 

Fig. 2. The central H2 number densities of Type 1 (OH94, blue 
solid line) and Type 2 (LD01, red solid line) critical MBESs as 
a function of core mass. The green solid line represents central 
densities of critical BESs corresponding to the mean tempera- 
tures of the Type 1 MBESs. Also plotted are the outer number 
densities of the respective cores types (dashed lines). 




Core mass [M mn ] 



Fig. 3. The temperature structures of Type 1 (blue) and Type 
2 (red) critical MBESs as a function of core mass. The upper 
bounds of the shaded areas represent the outer temperatures of 
the cores, the lower bounds represent central temperatures. The 
solid lines indicate mean temperatures (see text). 



2 core. The temperature difference between the centre and the 
outer boundary is clearly larger the Type 1 core (about 2.5 K) 
than for Type 2 (about 1 K). This is caused by the different dust 
models. 

3.2. The central densities and temperature profiles of critical 
non-isothermal cores 

Figure [2] plots the central number densities of critical MBESs 
and BESs as functions of the core mass. Also plotted are the 
MBES number densities at the outer boundary. The overall ten- 
dency seen in this Figure that the central density decreases 
steeply with the core mass can be understood from eqs. Q and 
(O, according to which the central density (represented by the 
parameter A) is approximately proportional to M~ 2 when the 
non-dimensional radius £ out does not change. The slight de- 
viations from the M~ 2 power-law are caused by temperature 
changes and small variations of £ out discussed below. The cen- 
tral densities of BESs are higher than those of MBESs in the 
entire mass range, and the Type 1 cores are always denser than 
the Type 2 cores. The ratio of central and outer number densi- 
ties of the MBESs remains approximately constant, regardless 
of core mass. There is however a small increase in this ratio to- 
ward small core masses; we will discuss this in Sect. 13.31 

Figure [3] plots the temperature structures of Type 1 (blue) 
and Type 2 (red) critical MBESs as a function of core mass. The 
upper bounds of the shaded areas represent the outer tempera- 
tures of the cores, the lower bounds represent central temper- 
atures. Also plotted are the mean temperatures of the cores as 
defined by eq. ( lB.4t (solid lines). The temperature difference be- 
tween the core centre and the edge increases towards smaller 
masses. This is caused by the increase in the central density (see 
Fig-EJ which leads to greater attenuation of external radiation. 
Low-mass MBESs have noticeable temperature gradients while 
cores with larger masses are almost isothermal (Type 1 cores still 
present a temperature gradient of ~ 1 K at M — 5 M ). Looking 



at the mean temperatures, we see that for high mass cores the 
mean temperature approaches the outer temperature, indicating 
that most of the core mass lies in the outer parts of the core. For 
small core masses, the mean temperature approaches the central 
temperature, and the medium is more centrally concentrated. 

3.3. The critical non-dimensional radius fi 

The non-dimensional outer radius £ out characterizes the shape of 
the density profile which can be compared with high-precisio n 
observations like in the famous case of B68 (lAlves et al.ll200l . 
The critical value of £ out where an MBES becomes unstable 
against an increase in the external pressure is denoted here by %\ 
(see Sect. 12.3b . Figure |4]plots the value of £1 calculated accord- 
ing to (eq.|9]l as a function of core mass. The blue and red lines 
correspond to Type 1 and Type 2 models, respectively. Also plot- 
ted in the Figure is the isothermal critical value £0 ~ 6.45 (green 
dashed line), which does not depend on core mass. 

One can notice here that for Type 1 cores %\ increases clearly 
towards smaller masses, whereas for Type 2 cores £1 is nearly 
constant and close to isothermal value £0, except for the wig- 
gle at the low-mass end. The different behavior of £1 in the two 
cases is due to the different temperature structures of the two 
core types (Fig.[3j. Consider first Type 1 cores. At high core 
masses, the temperature gradient of the critical core is ~ IK, 
whereas for small core masses the gradient is ~ 3.5 K. For Type 
2 cores, the high mass temperature gradient is ~ 0.3 K, whereas 
the low mass gradient is ~ 1.0K. Comparing with Fig. [4] we 
see that the critical radius correlates fairly well with the behav- 
ior of the temperature gradients; for Type 1 models, £1 grows 
as the temperature gradient grows, and for Type 2 models, £1 
is approximately constant until the temperature gradient starts 
to grow toward small core masses. A more general conclusion is 
that small mass cores with high central density and steep temper- 
ature gradient are able to remain stable for slightly higher values 
offi. 
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Fig. 4. The non-dimensional critical radius %\ for Type 1 and 
Type 2 models (blue and red lines, respectively) as a function of 
core mass. Also plotted is the isothermal critical radius £o (green 
dashed line). 



The increasing temperature gradient towards lower masses is 
also reflected in the density contrast, 



Pc 
Pout 



T, 



(10) 



For MBESs the critical density contrast is generally close that of 
a critical BES (~ 14), but exceeds the isothermal value for the 
smallest masses. 

We note that the values of £i plotted in Fig. [4] are subject 
to numerical uncertainty, mostly due to Monte Carlo fluctuation 
in the temperature calculations, which can be up to 0. 1 K. This 
fluctuation is reflected on the pressure derivative (eq.|9j through 
its dependence on /3, A and iff and their numerical derivatives. To 
analyze the effect of numerical fluctuation, we performed the it- 
erative calculations presented in Sectf2]with multiple £ out grids 
of varying accuracy (meaning specifically the spacing between 
neighboring values of £ ut)- It was found that the values of ^\ for 
a given core mass can change by up to 0. 1 depending on the grid 
used, but the general shapes of the lines plotted in Fig.|4]remain 
more or less the same in all cases studied. The numerical fluc- 
tuation is probably responsible for the small dip in the Type 2 
data present in Fig. [4] Furthermore, we would expect £i to ap- 
proach £o in the large core mass limit (nearly isothermal cores) - 
the apparent lack of convergence is here attributed to numerical 
un certainty. 

iGalli et aTl d2002h studied the structure and stability of 
MBESs by varying the external pressure and the strength of the 
ISRF (assuming outside shielding corresponding to Ay = 1). 
They report different values for the density contrast than the ones 
presented here. In the particular case of M — 5 M , they derive 
Pc/Pout ~ 20. This is an interesting result, as our models predict 
nearly isothermal values for the density contrast for high core 
masses. We attempted to reproduce their results by considering 
the density profile for the M — 5 M Q sphere as calculated here, 
but calculating the dust temperature assuming Ay = 1 . We cal- 
culate d the gas temperature in the same way as in IGalli et all 
d2002l) . i.e. balanc ing heating and c ooling functions using the 
parametrization of iGoldsmithl d200lh . 
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Fig. 5. The physical radii of Type 1 (blue solid line) and Type 
2 (red solid line) critical cores as a function of core mass. The 
black curve represents homogeneous spheres in virial equilib- 
rium (zero outside pressure) with temperature corresponding to 
the mean temperature of the 1 M Type 1 MBES. Also plotted 
are the physical radii of critical isothermal cores corresponding 
to the mean temperature of each non-isothermal core (dashed 
lines; see text). 

at the core center, decreasing to ~ 9 K at about half radius. At the 
core edge, the temperature rises again to above ~ 9.5 K. Using 
the gas temperature profile to solve the modified Lane-Emden 
equation (eq.|6]l and solving the density contrast from eq. dTOb 
yields p c /p out ~ 20. 

The high density contrast is probably due to the increased 
temperature at the core edge, where the increased pressur e sta- 
bilizes the core. The parametrization of lGoldsmithl d200lh does 
not take into account the geometry at the edge of the core; in 
a more "realistic" scenario, increased photon escape probability 
at the core edge may allow gas temperatures to keep decreasing 
toward the core edge unless the cloud is so weakly shielded that 
the external UV field is able to heat its surface layers (through 
photoelectric effect). T his scenario has recently been studied by 
lJuvela & Ysardl d201 lb . In the case of lower gas temperatures at 
the core edge, the predicted density contrast is probably closer 
to the isothermal value. However, proving this statement would 
certainly require a quantitative study. 

Finally we note that while there is dependence on core mass, 
our models predict that the values of %\ are equal to the isother- 
mal critical value £o ~ 6.45 up to an accuracy of ~ 5%. The 
density contrasts of the critical MBESs are also very similar to 
the density contrast of the critical BES, ~ 14. 

3.4. Physical radii and boundary pressures 

The physical radius, R out , of an MBES of given mass depends 
on the central temperature and the non-dimensional radius £ out 
according to 



^nnt — 



GM 



(id 



We calculate T„ 



12K 



where c s c = kT c /m is the sound speed in the core centre. For 
critical spheres, £ out = fi is roughly constant (Fig.|4j, and 
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Fig. 6. The boundary pressures p/k of critical Type 1 (blue line) 
and Type 2 (red line) MBESs as a function of core mass. Also 
plotted are the boundary pressures of BESs with temperatures 
corresponding to the mean temperatures of the MBESs (dashed 
lines). 



changes in the central temperatures are modest, so the radius 
is roughly proportional to the mass. 

We plot in Figure|5]the physical radii of critical Type 1 (blue 
solid line) and Type 2 (red solid line) cores as a function of core 
mass. The critical radii of the "corresponding" BESs with the 
temperature T — T, the average temperature of the MBES of the 
same mass, are plotted with dashed lines. The radii represent the 
minima of the static equilibrium below which compression leads 
to the gravitational collapse. Also plotted is the physical radius 
of a virialized homogeneous, isothermal sphere with zero outer 
pressure, for which R = GM/5c;r. Also in this case we have 
used the assumption T = t. One can see that the physical radii 
of MBESs are only slightly larger that those of the corresponding 
BESs. For the latter the critical radius can be approximated by 
R * GM/lAc;. 

To conclude this Section, we make a note about the bound- 
ary pressures of the critical MBESs. The boundary pressures 
p/k of Type 1 and Type 2 cores as a function of core mass are 
plotted in Fig. [6] (blue and red lines, respectively). The bound- 
ary pressures of the corresponding isothermal spheres are plot- 
ted with dashed lines. These diagrams represent the maximum 
pressures above which an increase in the outer pressure leads 
to a gravitational collapse. For sub-critical cores, i.e. those with 
fout < £ii the boundary pressure is lower and the physical ra- 
dius is larger. The boundary pressures of critical Type 2 cores 
are smaller than those of Type 1 cores. Low-mass critical BESs 
have higher boundary pressures than the corresponding MBESs, 
but the difference decreases towards larger masses and is practi- 
cally zero for M = 5 M . 

We have com pared our calculations with those of 
iKeto & Field! d2005l) . who presented in their Fig. 1 a p-V curve 
for an MBES of M - 5 M . For the thermally supported 
sphere. Ike to & Field! (120051) calculate for the critical sphere val- 
ues of the gas density and boundary pressure of ~ 10 4 3 cm~ 3 
and ~ 5.5xl0 3 cirT 3 K, respectively. As also indicated by Figs. [2] 



and|6] we predict ~ 10 41 cm -3 and ~ 5xl0 3 crrT 3 K for the cen- 
tral density and boundary p ressure, respectively. Our results are 
thus comparable to those of lKeto & Field! d2005l) . 

4. Discussion 

In this Section, we discuss the implications of the results pre- 
sented in the previous Section. We make some general remarks 
about the nature of the stability and discuss the implications of 
the physical properties of the critical cores on the core chem- 
istry. We conclude the Section with a discussion on very small, 
roughly Jupiter mass cores. 

4.1. The nature of the stability 

iBonnor! (Il956l) showed, by studying the second pressure deriva- 
tive, that the minimum volume for a stable core corresponds to 
the configuration for which the pressure derivative with respect 
to volume dp/dV at the core boundary becomes zero. The result 
applies here as well, although we consider non-isothermal cores. 
It follows (see also Sect. l3.4b that each physical radius shown in 
Fig. (|5]i represents the smallest stable core configuration for a 
given core mass and model type. That is, for a given core mass, 
configurations with R > Ri are stable, whereas all configurations 
with R < Ri are unstable. 

At this point, however, it is prudent to consider what the sta- 
bility actually means. The stability condition (eq.[9j> deals with 
small first-order fluctuations of the relevant parameters (central 
density, central temperature and non-dimensional radius). Thus 
in this context, the statement that a core is stable implies that 
the core is stable against small fluctuations of these parameters. 
The possibility that sufficiently large fluctuations could induce 
collapse even in a "stable" core cannot be ruled out based on the 
analysis presented here. 

Also, our analysis does not give insight into what happens af- 
ter a core begins to collapse. It may well be (as has been demon- 
strated in the literature and by our preliminary MHD collapse 
models) that "unstable" core configurations settle into a stable 
state relatively soon after collapse begins. This is due to the in- 
crease of the temperature inside the core: as the core collapses, 
the thermal pressure inside the core increases and may be able to 
support the slightly contracted core, depending on the core cool- 
ing mechanisms. Indeed, the stability condition presented here 
provides an initial condition for the collapse of non-isothermal 
cores, rather than trying to predict the "absolute" stability of the 
cores. 

4.2. Physical structures 

Figures Q] to [3] demonstrate that a Type 1 MBES has a different 
physical structure than a Type 2 MBES of the same mass. In gen- 
eral, Type 1 cores are denser, warmer and smaller (Pig.[5j than 
Type 2 cores, i.e. Type 1 cores are more compact. Toward small 
core masses, the differences in radius and density diminish, but 
Type 1 cores always present a larger temperature gradient. 

If one computes the mean temperature f of a given MBES, 
a "corresponding" BES of the same mass can be constructed. 
Figure [2] shows that the central density of a BES constructed 
in this way is larger than that of the MBES - this fact is also 
apparent in Fig.[TJin the special case of M — 0.25 M B . The BES 
is, however, smaller. This result is general for small core masses. 

2 The central number density plotted in Fig.[2]is the H 2 density, so 
that the fofai hydrogen number density is twice that plotted in Fig.|2] 
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However, when using the gas temperatur e the situation chang es 
to the opposite for larger masses (see e.g. lKeto&Fiei3l2005l) . 

One can, of course, construct a BES with outer radius R ex- 
actly equal to that of a given MBES. In this case, the tempera- 
ture of the BES is clearly lower than the mean temperature of the 
MBES, but the density structures turn out to be nearly identical, 
with very small differences caused by the different temperature 
structures. 

4.3. Chemistry 

We now briefly discuss the impact of the physical properties of 
the cores on their chemical evolution. 

Let us first concentrate on the differences between the Type 
1 and the Type 2 MBES. As discussed above, a Type 2 MBES 
is larger than a Type 1 MBES of the same mass, although the 
former is less dense and colder. The typical density difference in 
the center is however less than 15% between the core types, so 
in this case, the major factor in the possible different chemical 
evolution of the cores is the temperature gradient. 

Even small temperature changes can have significant con- 
seque nces for the chemis try and the interpretation of observa- 
tions. [Ajkawaetal ] (l2005h studied the chemical evolution in col- 
lapsing clouds with initial conditions close to critical Bonnor- 
Ebert spheres. Comparing models with a central density of 
3 x 10 6 cirT 3 and temperatures in the range 10K-15K, the 
molecular depletion in the central parts of the core was found 
to be very sensitive to the temperature. In particular, the abun- 
dances of NH3 and N2H + could differ by an order of magnitude 
or more. However, because the central abundances are always 
low compared to the outer regions, it may be impossible to ob- 
servationally confirm such a temperature dependence. 

Let us then turn our attention to the differences between the 
MBES and the BES. As discussed in Sect.gH the BES is al- 
ways more compact than the corresponding MBES. The temper- 
ature dependence is more complex than when comparing two 
MBESs - the BES is warmer in the center, but colder toward the 
edge. These differences may again produce notable differences 
in chemical evolution. 

It is difficult to determine core temperature profiles from 
observations and, similarly, the oretical pre d iction s contain sig- 
nificant sources of uncertainty. iGoldsmithl d2001l) showed that 
molecular depletion can significantly increase the gas tempera- 
ture predicted for cores. The effect can be several degrees even 
in dense clouds although, once the density is ~ 10 5 ctrT 3 or 
above, r gas tends towards the du st temperature irrespect ive of the 
decreasing line cooling power. lJuvela & Ysardl J20Tlt) demon- 
strated further how, in addition to abundance variations, the ra- 
dial temperature profiles can be modified by the velocity field 
(affecting the efficiency of line cooling) and the grain size distri- 
bution (affecting the gas-dust coupling). Therefore, each hydro- 
static object is likely to exhibit slightly different density profiles. 

The above implies that chemistry plays a role also in the 
dynamical evolution of a core. A study of the chemistry of the 
MBES including line cooling is planned - this will allow a more 
realistic estimate of the gas temperature, also serving as a start- 
ing point for dynamical studies. In this context, one can also 
quantitatively study how the gas temperature affects the stability 
analysis presented here. Finally we note that the stability analy- 
sis presented in this paper indicates that an MBES can be just as 
stable as a BES of the same mass and (roughly) the same size, 
making it possible for advanced chemistries to develop in both 
cases in a typical core lifetime. 



4.4. Jupiter mass non-isothermal cores 

In the above, we present results of a stability analysis for core 
masses ranging from 0.1 M to 5.0 M e . It has been suggested 
that small dark clouds of roughly Jupiter mass could make up 
for part of the submm sources observed in e. g. SCUBA maps , 
or even account for part of the dark matter ( Lawrence] 1200 lb . 
Alt hough this scheme has later been found to be unlikely (see 
e.g. iDrake & Cooll2003t lAlmaini et al1l2005l) . the existence of 
such objects is not entirely ruled out. For this reason, we thought 
it could be interesting to model these roughly Jupiter mass 
clouds (hereafter JMC) as either a BES or an MBES, to find 
out what kind of constraints would be imposed on the stability 
of these objects. 

In our calculations thus far, we have assumed Ay — 10 at the 
core edge, so that the objects are embedded in a bigger construct, 
e.g. a molecular cloud. However, if the JMCs are thought to be 
isolated, then the assumption of Ay = 10 at the core edge is 
no longer justified. Indeed, if a JMC is isolated, the outermost 
parts are most likely ionized and heated to high temperatures 
by the ISRF, and molecular gas can be thought to be found at 
a minimum visual extinction of Ay ~ 2. We have thus carried 
out the stability analysis for a core mass of 0.01 M (roughly ten 
Jupiter masses), assuming outside Ay = 2. 

As a point of reference, let us first look at the BES. The crit- 
ical radius of a 10M Jup BES at 7 K is 140 AU, and the corre- 
sponding boundary pressure is p/k = 4 x 10 8 ctrT 3 K. As dis- 
cussed in Sect.|3] these values represent the minimum physical 
radius and maximum boundary pressure for a core in stable equi- 
librium. The critical radius is directly proportional to the mass 
and the critical boundary pressure is proportional M~ 2 . So for 
example, a 1 Mj up BES at 7 K would have a critical radius of 
1 4 AU and a bound ary pressure of 4 x 10 10 cm -3 K. According 
to|L awrenc the blank-field submm sources have angular 

radii of about 1" or less, and if interpreted as cold dark clouds 
they should be nearby objects with a characteristic distance of 
100 pc. Hence a typical radius for this kind of cloud should be 
about 100 AU. This radius roughly equals the critical radius of a 
10 M Jup BES. For a 1 M Jup BES at 7K, the radius R = 100AU 
corresponds to a sub-critical dimensionless radius of £ out ~ 1. 
The boundary pressure for this configuration is still very high, 
p/k~2x 10 8 ctrT 3 K. 

Let us then look at the MBES. For a critical Type 1 core, the 
physical radius is ~ 200 AU. The temperature gradient is rather 
steep (from 5K at the center to 10.5 K at the edge), and the 
boundary pressure is p/k ~ 10 8 ctrT 3 K. The temperature gra- 
dient of the Type 2 core is shallower (from ~ 4.8 K at the centre 
to ~8.5K at the edge), but R and p/k are nearly identical to 
those of the Type 1 core. This analysis suggests that the mini- 
mum size of a stable 10Mj up MBES is ~ 200 AU, irrespective of 
the dust model used. However, the boundary pressure required to 
maintain the critical configuration in equilibrium is unlikely to 
be present in interstellar conditions, at least if the core is isolated. 
For clearly subcritical configurations (£ < £1, R ~ 1000 AU), the 
boundary pressure is still > 10 7 cm" 3 K. 

We conclude that due to the high boundary pressures re- 
quired, roughly Jupiter mass critical BESs or MBESs are un- 
likely to be able to exist as isolated objects in interstellar 
space where the pressure is typically ~ 10 4 - 10 5 ctrT 3 K (e.g. 
iMcKee & Ostrikerl 120071) . Subcritical configurations can have 
radii of ~ 1000 AU, but the boundary pressure required for stable 
equilibrium is very high also for these. 
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5. Conclusions 

We studied the stability of non-isothermal Bonnor-Ebert 
spheres. The physical parameters of critical cores were derived 
for a range of core masses. Two differe nt types of dust grains, 
corres pondi ng to optic a l dat a from lOssenkopf & Hennind 
(fl99l and iLi & Draind (EooTh . were considered in the mod- 
elling. The analysis was accordingly separated into two distinct 
core types (Type 1 and Type 2, respectively), each containing 
one type of dust grains. 

As a general trend, the central density of a core increases 
as core mass decreases. The increase in density also increases 
the temperature gradient, because of more efficient attenuation 
of the external radiation field (in this paper, we study the dust 
temperature only). However, the absolute values of the temper- 
ature depend on the adopted dust model. Type 1 cores present 
larger temperature gradients than Type 2 cores, but Type 2 cores 
are colder. The difference also translates to core size: due to 
the larger temperature, the thermal pressure in the centres of 
Type 1 cores is larger and hence a Type 1 core can assume a 
smaller size than a Type 2 core of the same mass. Thus, Type 
1 cores are more compact. Isothermal cores corresponding to 
non-isothermal cores (i.e. with same mass and same average 
temperature) are more centrally dense and slightly smaller than 
their non-isothermal counterparts. Considering gas temperature 
instead of dust tem perature may affect core stability (see e.g. 
lKeto& Field! l200l . Our results should, however, hold well at 
least for cores with masses up to ~ 0.5 M , for which the aver- 
age density is ~ 10 5 cm 4 . 

It was found that the non-isothermal critical radius %\ in- 
creases above the isothermal critical value £0 ~ 6.45 to- 
ward small core masses. The effect is particularly pronounced 
for Type 1 cores, which present large temperature gradients. 
Nevertheless, the change in i;\ with core mass is rather small, 
and in the mass range studied here, %\ is equal to £0 to within 
5%. Furthermore, the ratio of central and outer densities for crit- 
ical MBESs was found to present a similar increase toward small 
core masses, but the ratio is also in this case close to the isother- 
mal value (~ 14). 

Although the differences in physical parameters between the 
two MBES types are not great, the different temperature gradi- 
ents in particular may affect chemical evolution in these objects, 
which might be observable through line emission radiation from 
these objects. This may also the case when comparing isother- 
mal and non-isothermal cores. Line radiation cooling can also 
affect the stability of the cores. A quantitative study of chemical 
evolution in the different types of cores could thus be justified. 

We studied the physical parameters of MBESs of roughly ten 
Jupiter masses (~ 0.001 - 0.01 M ) in order to validate or dis- 
qualify their possible existence as isolated objects. It was found 
that the boundary pressure required to maintain critical cores 
equilibrium is unlikely to be found in interstellar conditions. 
Even subcritical configurations require boundary pressures of 
p/k ~ 10 7 ~ 8 cirT 3 K. We conclude that very low mass MBESs 
are unlikely to be able to exist in the interstellar medium as iso- 
lated objects - this applies to very low mass BESs as well. 

Finally, we note that the stability analysis presented in this 
paper considers small, first order variations of the relevant pa- 
rameters (central density and temperature, core radius). Thus, 
stable cores as defined by the stability condition derived here 
are stable against linear perturbations, and sufficiently large per- 
turbations could induce collapse even in these "stable" cores. 
Furthermore, what happens after collapse begins is not predicted 
by the equations. However, the near-critical gas spheres studied 



here represent plausible, albeit idealistic models for cores at the 
very beginning of collapse. 
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Appendix A: The stability condition 

In this Appendix we present a short derivation of the general 
form of the derivative of the boundary pressure with respect to 
core volume (eq.|9j. We also present an expression for the total 
number of particles in terms of f3, A and £ out . Throughout the 
calculations, we make use of equations (f2]i to (O in the main 
text. 



A. 1. The pressure derivative 

Calculating the variation of the pressure (eq.|4j) with respect to 
the three non-dimensional variables /3, A and £ out yields 



dp = -5(pT) = A t 7iG5(AI3e-' lJ ) 



= AnGpAe^ 



A- l 6A+p- l 6/3- 



JdijA 



The variation of the cloud volume is 



5V = Anr^dr 



2^ Vout« 



1 _! 

-A t; oyA 5A + <5£ out 



(A.l) 



■(A.2) 
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The above equations then yield the pressure derivative in general 
form: 



dp 

— = 2G/T 
6V H 



1/2,5/2^-3 



. (A3) 



Identifying 2G/T I/2 ,l 5/2 £ ra 3 t e _ ' A = 2p/3V leads finally to eq. ©. 



-ATlnZ, and the entropy can be obtained from 
Below, in the derivation of the internal energy, 



is then F = 

s = -(f) 

\dT Jn,V 

we use a more straighforward method. In any case, Q has to be 
determined from the density distribution, and in this prescription 
it is a function of /V (mass), V (outer radius), and an additional 
variable £ out which describes the distribution of mass within the 
volume V. 



A.2. Mass conservation 



B. 1 . Internal energy U 



We now present a non-isoth ermal version o f the mass conser- 
vation equation discussed bv lBonno r! (ll956l) . Using eqs. (f2]i and 
(0, we can write the total number of particles as 



4tt r 

N=— pfdy 
m Jo 



An 



3/2 -,-1/2 



-fA 



(A.4) 



The mass of the core is M = Nm. Since N is calculated for a 
single core (that is, for a given value of £ utX one can move /3 
and A out of the integral; for any single core, these are constant. 
In the integral, the variable £ represents the internal structure of 
the core and hence we make use of the solution to the Lane- 
Emden equation (eq.[6j> that is unique to each core. The above 
formula simplifies to 



N=-/3 i/2 A- 1 ' 2 f To 
m 



dip 



(A.5) 



where r out represents the temperature contrast of the core, that is 
the ratio of outer and central temperatures. During the collapse 
process, we expect the total number of particles (i.e. the total 
mass of the core) to be conserved, i.e. that 6N = 0. Calculation 
of this condition leads to the equation 



A- l SA=p- l 6p + ^e-^) 



(A.6) 



It is readily verified that in the isothermal limit (5/3 = 0, t = 1), 
eq. dA.6b reduces to its isothermal analogue (eq. 2.11 in lBonnorl 
1956). 



Appendix B: Thermodynamics of the MBES 

iBonnorl d 19561) discussed the form of Boyle's Law for a mas- 
sive isothermal gas sphere. The radial change of the temperature 
causes a slight modification to Bonnor's equation of state, and 
we think there is a reason to make a brief excursion to the ther- 
modynamics of an MBES. Besides writing down the equation of 
state, we also present the first law of thermodynamics in the case 
of an MBES which can be useful for d ynamical calculation s or 
when discussing thermal instability (see lKeto & Fieldll2005l) . 

It turns out that the internal energy of an MBES can be de- 
scribed by the simple ideal gas expression of internal energy 
when the temperature T is replaced by the mass-averaged tem- 
perature f defined below. In general terms, an MBES can be 
considered as a thermodynamic system containing N molecules 
in a volume V, at the temperature T . The thermodynamic po- 
tentials and the equation of state can be derived much in the 
same way as for an ideal gas (see e.g. lLandau & Lif shitzll 1 969t 
lMandll[l988h from a partition function Z which is the product 
of the perfect gas partition function, Z p (N, V, T) (which may 
also contain rotational and vibrational parts), and the configu- 
rational partition Zn = exp(-Q./ kT), where Q. is the gravita- 
tional potential energy. For example, the Helmholtz free energy 



Assuming that the medium consists of ideal gas with effec- 
tive degrees of freedorrQ / = 3, we c an equate internal en- 
ergy with total thermal ener gy (see e.g. IChandrasekhari [19571 : 
iKippenhahn & Wei gertlll 99*41) and write 



U 



= 6n f 

Jo 



pr 2 dr . 



(B.l) 



Using eqs. to ©, we can write the integral in terms of f3, A, 
and f , and obtain 



/"foul 

U = 24n 2 G /3 5/2 A- l/2 J £ W<f . 



(B.2) 



On the other hand, the assumption of ideal gas means that the 
total internal energy of an isothermal core is 



U = -NkT . 

2 



(B.3) 



It can be shown that eq. ( IB. 31 ) is valid also for a non-isothermal 
sphere if T is replaced with the average temperature t defined 



T 5 



Io T PdV 

JoPdV £out T out(f) 



(B.4) 



When substituting the expressions of T and N from eqs. ( IB. 41 ) 
and iAAi to eq. ( IB. 31 ) one obtains the formula for U given in 
eq.(lB~2b. 

B.2. Gravitational potential energy 0. and the virial theorem 

The total potential energy of a gas sphere of radius R can be 
written as 



£2 



C gm j*, a r GM p 3, 

dM = -An — T-rdr, 

Jo r Jo r- 



(B.5) 



where p = p(r) and we have written the mass of a spherical shell 
as dM = Anr 2 pdr. Using the hydrostatic equilibrium condition 



dp _ GMp 
dr r 2 



(B.6) 



we can rewrite eq. ( IB. 51 ) into an easily integrable form: 

Q = 47rf — r 3 dr = 4npR 3 - I2n [ pr 2 dr. (B.7) 
Jo dr Jo 

The latter term in eq. ( IB. 7b can be identified as twice the thermal 
energy (see eq. lB.U . Rearranging terms yields the virial theorem 



(B.8) 



2U + £l = 4npR 3 = 3pV . 

3 The very cold (T ~ 10 K) ambient temperature in prestellar cores 
effectively removes the rotational degrees of freedom, so that even H 2 
can be treated as a monatomic gas, at least to reasonable accuracy. 
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Fig.B.l. The function /(&ut) for a BES and an MBES of M = 
IM Q . 

Writing Q. in terms of A and £ out and using eqs. OJ and ®, it 
can be shown that 



-^/(W = -(|) 1/3 GM^-^ 0Ut)! 



(B.9) 



where / is a function of tout only. To illustrate, we plot in 
Fig. ( tBTTb the ratio /(&ut) = Q/(-GM 2 /#) for a BES and an 
MBES of M = 1 M . Also marked is the isothermal critical 
value to ~ 6.45. The Figure shows that the gravitational po- 
tential energies of the BES and the MBES differ from each other 
somewhat. For small values of £ out the cloud is nearly homoge- 
nous and /(^out) approaches the value 0.6 = 3/5. 

6.3. Equation of state 

Substituting the expression of the internal energy in terms of the 
average temperature T to the virial theorem (eq. IB.8l ) one obtains 
the equation of state at the core boundary 



- Q 
pV = NkT + y . 



(B.10) 



Using eq . ( 1B.9I >. this can be written in the form suggested by 
iTerletskvl (Il952i Bonnor's eq. 1 .2), criticized by Bonnor, if the 
numeric factor a is taken to be a function of tout- 



The validity of this formula can be verified by using on the right 
hand side one of the Maxwell relations 



dS_\ _ ldp\ _ Nk 



dV 



dT) 



V.N 



V 



On the other hand, from the definition of heat capacity: 



-C =- — 

T V 2 T 



(B.14) 



(B.15) 



where we used the result that for a monatomic ideal gas 
C v - \Nk. Equations ( IB. 14b and ( IB. 15b together yield 



NkT 3 
TSS = 6V+-NkST. 

V 2 



(B.16) 



Substituting this to the right hand side of eq. ( IB. 131 ) one obtains 
the identity 



6U = -Nk ST , 
2 



(B.17) 



in accordance with eq. ( IB. 31 ). 

The total energy of the core equals to the sum of internal en- 
ergy and gravitational potential energy, E — U + CI, from which 
follows that 



5E = 6U + SO. . 



(B.18) 



According to eq. (IB.9I ) we can write O = Q(V, £ out ) when the 
mass is kept constant. Then 



dQ. dQ. 
SQ = — SV + —dt 0U t ■ 
oV <9£,ut 

Plugging eqs. (IB. Ill ) and (IB. 191 ) to eq. ( IB.18I ) yields 
dQ. 

SE = TSS - pSV + — Sto* . 



(B.19) 



(B.20) 



Using eq. ( IB.16I ) above the change of the total energy takes the 
form 



3 - Q d£l 
SE = -NkST - —SV + — S^ut ■ 
2 3 V dtout 



(B.21) 



The above equation is another expression of the first law of ther- 
modynamics. 



B.4. First law of thermodynamics 

Considering the internal energy only, the first law of thermody- 
namics states that 

SU = TSS - pSV - ^SV , (B.l 1) 

oV 

where the last term represents the work done by the gravitational 
potential energy. For a reversible process, the first term can re- 
placed by heat dQ supplied to the system. From eq. dB.91 ) it can 
be seen that 

dV 3V v 

Furthermore, using the equation of state the change of the inter- 
nal energy becomes 

SU=TSS-?^SV. (B.13) 



